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ENGINEERING DESIGN HANDBOOK 
TABLES OF ‘THE CUMULATLVE LTNOMIAL PROBABLLITIES 


FOREWORD 


The Binomial Distribution and the Incomplete Beta-Function Ratio (or Pearson Type I 
distribution) are probability distributions which are of great importance in both the theory of 
probability and in its applications. 


The binomial probabilities which are generated from (q +P)”, where n is an integer and q tp = 1, 
constitute a discrete probability distridution which has very wide application. The Incomplete 
Beta-Function Ratio, defined by 


XO P Beas ety Blan oF 1) 
f po py { pod py) dps erten ctl) 
0 “0 ; Blown - eo) : 


which has been tabulated by Pearson’ is continuous. Nevertheless, precise values of one distribution 
can be derived directly from the other by means of the simple identities 


Pump) = S ("\or q" vs ilen etl) 


si o 
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Se a. ee . 
~P(ct imp) = 2 ("yp qe rsd, pn- cet h 
7-6 


=1- {ke t+) c¢) 


Since the transition from the binomial distribution to the Incomplete Beta-Function requires 
only a simple change in parameters, which function one should tabulate is theoretically 
unimportant. The choice, however, is one of considerable practical importance. The Tables of the 
Incomplete Beta-Function' are so arranged that in many applied problems a preat deal of page 
turning is required to extract the needed probabilities. It is believed that the tabulation of the 
binomial distribution in the format of the following tables is favored by many applied statisticians, 
especially in view of the widespread use of binomial probabilities in industrial and Government 
problems, in sampling inspection, and in statistical analyses generally. Indeed, the format used 
herein will in many practical problems effect a saving in time and provide some convenience in 
searching for the required probabilities. 


These tables arose from computations related to the solution of some military problems during 
the late 1940's and early 1950's ai the Ballistic Research Laboratories of sche Ordnance Corps. U.S 
Aimy, Aberdeen Proving Ground, Maryland. Only the cumulative binomial probability for single 
*This pamphiet replaces AMCP 706-109. Tebies of the Cumulative Binomial Probabilities, Vanuaty 1971. 


' Superscript numbers refer to References listed 
on p. xxxii. 
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chances of 0.5 or less is tabulated. The Introduction points out transitions for obtaining the 
probability values for the upper half of the range (0.5 to 1.0) from the tabulated values of the lower 
half. ~ 


The computation of the tables was originally suggested by and carried out in collaboration with 
Dr. Frank E. Grubbs; the tabulation was done on IBM Relay Calculators at the Ballistic Research 
Laboratorics during 1952. Valuable contributions were made to the work by Dr. Max Lotkin and 
others of the Computing Laboratory of the Ballistic Research Laboratories. In order to make the 
tables generally available to the scientific community, arrangements were made for the publication 
of the tables in an economical way through the Government Printing Office in 1953. 


In the past year, copies of the tables have become unavailable. For this reason, the tables were 
made a part of the Engineering Design Handbook Series, ana additional copies (with corrections) 
were printed. Hence, it appeared desirable to prepare a new Introduction that would be morc 
up-to-date in applications. The work on the draft of the new Introduction was borne largely by 
Professor William C. Guenther of the University of Wyoming. It is hoped that the new printing, 
therefore, will serve wider purposes. Those contributing to this revision are: Frank E. Grubbs, 
William C. Guenther, Leslie E. Simon, and Helen J. Coon. 


The Engineering Design Handbooks fall into two basic categories, those approved for release and 
sale, and those classified for security reasons. The Army Materiel Command policy is to releasc 
these Enginecring Design Handbcoks to other DOD activities and their contractors and other 
Government agencies in accordance with current Army Regulation 70-31, dated 9 September 1966. 
It will be noted that the majority of these Handbooks can be obtained from the National Technical 
luformation Service (NTIS). Procedures for acquiring these Handbeoks follow: 


a. Activities within AMC, DOD agencies, and Government agencies other than DOD having need 
for the Handbooks should direct their request on an official form to: 


Commanding Officer 

Letterkenny Army Depot 

ATTN: AMXLE-ATD 
Chainbersburg, Pennsylvania 17201 


b. Contractors and universities must forward their requests to: 
National Technical Information Service 


Department of Commerce 
Springfield, Virginia 22151 


(Requests for classified documents must be sent, with appropriate “Need to Know” justification, to 


Letterkenny Army Depoi.) 
Comments and suggestions on this Handbook are welcome and should be addressed to: 
U.S. Army Materiel Command 


ATTN: AMCRD-TV 
Washington, D.C. 20315. 
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1. INTRODUCTION 


The binomial probability distribution is of considerable importance in both the theory of 
probability and its applications. It is the appropriate model to use in computing probabilities fora 
series of events or trials, or for any experiment which has the following properties: 


(a) The result of each trial car be classified ito one of ty,o categories, e.g., success 
or failure, “good” or “defective”, conforming or nonconforiuing. 


(b) The probability p of a sv ccess is the same for each trial. (1-1) 
(c) Each trial is independent of all the others in the probability sense. 
(d) The number of trials is fixed at, say, 7. 


Such a series of events is generated by tossing an unbiased coin 5 times and regarding the 
occurrence of a head as a success (p = 1/2), or by rolling an unbiased die 10 times and regarding the 
occurrence of a six as a success (p = 1/6), etc. In either case, it is easy to verify that Conditions I-1 
are satistied. In a number of situations in which one desires to use the binomial distribution, there 
may be some doubt that conditions (b) and (c) are really met. Frequently in such cases, however, 
the experimenter can control certain factors in the performance of the series of experiments which 
make the Binomial Conditions 1-1 realistic in practical situations (as illustrated in Examples 1-1, 

-2, and 1-3), 


In the die example, the same die should be rolled each time. Exactly the same procedure or type 
of trial is carried out ten times. Each roll and, in general, each repetition is referred to as a “trial”’ of 
the experiment. The same terminology is used if ten dice (all unbiased) are thrown simultaneously. 
That is, cach die rolled is called a trial of the experiment even though ten different (but equivalent) 
“trials” have been performed. Thus, each trial is either a separate occurrence or a repetition in the 
same experiment. The kind of trials with which we are concerned here should satisfy Conditions 
1-1, 


Suppose that an experiment consisting of 7: trials satisfying Conditions 1-1 is performed and Ict X 
be the (random) number of successes obtained in the 7 trials. In most textbooks on probability, it is 
shown that the probability that X =x is 


? 
b(x np) =(") p* (1 —p)"—>*,x=0,1,2,...,2 (1-2) 
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‘ 
: The name binomial arises from the fact that dCx,1,p) is a term in the expansion of the binomial 
y (q +p)" or (Ci -p) tpl” where q= 1 — p. In most applications we necd a sum of terms rather than 
{ » aos . oY . . . 
: an individual term given by Eq. §-2. The accompanying table gives the right-hand sum 
ae ; 
Rd es 
: P(ea,p) = S bQnp) 
. i vte 
ane | 
a = Pr(X > 0c) (1-3) 
PS | = Pr(c or more successes inn trials) 


to seven decimal places for a = 1(1)150, p = .OLC.G1).S0 and also for p = .00i1¢.001).010. 


wii aia 
ad 


i? 


’ 
meat! errno. 


i 


We also define the left-hand sum 


Tae 





! 
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B(c,n,p) bQ,n,p) 
E 2 z (1-4) 


N 


PRX <c) 


Since the sum of the individual probabilities in a probability distribution must sum to unity, we 








have 
Bianp) = 1-— Me + Iynp) (1-5) iS 
The individual terms given by Eq. 1-2 are obtaincd by subtraction of successive tabled values. Thus, a 
a blenp) = Plemp) — Ple+ Vrep) (1-6) be 
a " It is unnecessary to tabulate sums for p > .50. When a p > .50 is encountered, one may ie 
aan interchange the role of p and | — p, ie., success and failure, and the desired sum may be obtained a 
: trom the accompanying tables for p < .50. The result is intuitively obvious but may be proved Ke 
algebraically. In Eq. 1-3 let x =n - y, which leads to & 
n-c ad 
if Plenp)= Si  b(y.n.1 ~ p)= Bln — e.n,1 = p) (1-7) i 
3 ye fe) 
oe Similarly, Eq. 1-4 yields 4 
Ay hi 
a re 
Ble.n,p) = Pn — ent — p) (1-8) k 
i 4 
a Some examples on the use of the tables follow. i 


Example 1-1 


In a 25-question, 5-answer multiple-choice examination, waat is the probability of getting 6 or 
more correct answers by (randomly) guessing? Exactly 6 correct answers? No more than 8 correct 


answers? 
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Solutions 
We observe that 


(a) Assuming that all the questions are answered, then each question answered is cither right or 
wrong, 


(b) For each question the probability of a correct guess is p = 1/5 = .20, 


(c) When each answer is a guess, presumably the answer to any one question does not influence 
the answer to any other question. (In many multiple choice type examinations there may be 
real doubt as to whether the independence condition is satisfied. Undoubtedly, an 
appropriate examination should be planned and constructed so as to make this condition 
scem as realistic as possible.) 


(d) The examination contains 25 questions. 


Thus, the binomial probability distribution is appropriate for application here with n = 25 and p 
= 20. Letting X be the number of correct answers, we find from the fables 


Pr ¥ > 6) = P(6,25,.20) = 3833106 
Pr(X = 6) = Pr(X 2 6) — Pr(X > 7) = P(6,25,.20) — P(7,25,.20) 
= 3833106 — .2199647 = 1633459 
Pr(X <8) = 1 ~ P(X 29) = 1 — P(9,25,.20) = 1 ~ 0467742 = 9532258 


Example 1-2 


Suppose that a machine produces items, 8 percent of which are defective. If a random sampte of 
40 items is inspected, what is the probability that 5 or more of the 40 items are defective? 


Solution 

Checking the Binomial Conditions 1-1, we observe that 

(a) Each item can be classified as a defective or nondefective. 

(b) Provided the machine remains in good condition, it is reasonable to assume that p, the 
probability that an item is defective, remains constant from item to item, at least for some 
suitably short interval of production. 

(c) If the machine is inclined to produce several poor items in succession, then the independence 
condition will not be satisfied. Therefore, in an attempt to make both (b) and (c) realistic, 
one might select the sample “at random” from some accumulation of items produced by the 


machine. 


(d) The experiment consists of “trials” of inspecting one item at a time until 40 items have been 
inspected. 
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é 


: . Thus, the binomial distribution seems to be appropriate with 2 = 40 and p = 08. Letting X be the 
‘ number of defectives observed in the sample, we find from the tables 


Pr{X > $) = P(S,40,.08) = 2132078 


- 


Example 1-3 
i For a certain variety of seed, it is hypothesized that 95 percent will germinate. If 50 seeds are 
a planted. what is the probability that 45 or less germinate? 
= | Solution 
“y 
Checking the Binomial Conditions 1-1 we observe that 
_ ( Asced either germinates or it docs not, 
ws fT 
{ (bv) If the seeds are selected at random from those seeds of the given type, it is not unreasonable 
‘ to assume that p, the probability that a seed germinates, is .95 for each seed in the sample. 


(c) It is reasonable to assume that seeds germinate independently of one another. 


t 
AE. 





(d) The number of ivials for the experiment is 1 = 50. 


Thus, the binomial distribution seems to be appropriate with 2 = 50 and p = 95. Letting X be the 
number of seeds that actually germinate in a sample of size 50, we want P(X < 45) = B(45,50,,.95) 
= 1 — (46,50,.95), Since p = .95 is not an entry in the table, we consider Y, the number of seeds in 
the sample that do not germinate. Now p = .0S and X¥ < 45 corresponds to Y = 5. Hence, the 


desired probability is 


ei are ee eT 





vai 
rm. 


PRY > 5) = P(5,S0,.05) = 1036168 
It is noted that substitution in Eq. 1-8 also yiclds the correct result. 


The following table might provide a useful summary: 
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2, RELATIONSHIP BETWEEN BINOMIAL PROBABILITY SUMS AND THE INCOMPLETE 
BETA-FUNCTION 


The Incomplece Beta-Function Ratio, 1 (A.B). Pearson,’ Harter. is defined by 
(ine = ae Pau O<u<! AB>O aI 
A,B) = , BAB) u (I — ua) du u A, (2-1) 


where (A.B) is the complete beta function of A and B, and 6(A,8) = MCA CBY/P U8) By 
repeated integration by parts itis easily shown that if 4 and B are integers, then 


1, (A.B) = PAA +B = yp) (22) 
Equivalently, we have 
Ion - +1) = Pleap) (2-3) 


Which is the probability of ¢ or more successes in a binomial type experiment, Moreover, the chance 
of ¢ or fess successes in 7 trials is given by 


I, -pQr--act hel - Ile+ la-aFl- Met lap) (2-4) 


In Table Lof the Tables of the Incomplete Beta-Function, Pearson’ has tabulated 1, (A.B) to 
seven decimal places for all combinations of p = .010.0191,00, 4 = SCS)11C1)50, # = SC.3)1 104950 
and A = B. (Some additional values are given in his Table HW) and Table IV.) Thus, for. and B 
integers every value tabulated by Pearson is available also in this binomial probability table. Since A 
+ B - J has a maximum value of 99 in Pearson’s table, then Eq. 2-3 can be evaluated only with 2’s 


up to 99 using his tables. 


Sometimes we desire a percentage point for a binomial distribution, i. ¢., given c, nm and a, we 
sack that p which satisfies either 


Plemp) = & (2-5) 
or 
Blenp) = a (2-6) 


In a number of practical problems linear interpolation in the binomial probability table for a 
percentage point will provide sufficient accuracy. If a higher degree of accuracy is desired, one 
could use any one of several tables. Leone, Haynam, Chu, and Topp’ have tabulated to eight 
decimal places the p which satisfies Eq. 2-6 for all combinations of a = 10(5)100, and aw = .0025, 
.005, OF, .025, .05, 1, .9, .95, .975, .995, .9975. In addition, one could use one of several tables of 
percentage points for the Beta Distribution. For example, Harter? gives to seven decimal places 


p = X(P,A,B) (2-7) 
for A, B = 1(1)40 and for twenty-three vasues of P or &. Here p is the solution of the equation 


1,(A,B) = P (2-8) 
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and consequently the solution of Ea. 2S with P= ae = dA anda oo +1 = B. The 2 ore entries for 
lharter’s table are 2 = .OOOL, .OUOS, QO! 005. .G1. 025, 08, 10 ¢.100.50. WP > 68. the one uses 
NAB) = ONCE PB UA). Pearson and Hardey.” Vogler and Thompson® have also prepared 
tables of the percentage points of the incomplete beta distributions. The only entries that these 
tables contain which are net used by Harter are 4 = 60, 8 = OO. P= Sand P= .75. 


As is well known, the percentage points of the binomial distiibution, or of the Incomplete 
Beta-Funetion Ratio, are equivalent to confidence limits on the true success ratio py having observed 
¢ successes ina fixed sample of size nv. For example, the value of p. equal to Py, aay, for which 


Bley, ) = hn, (e tla elhze 


gives an upper fa confidence limit on the true (unknown) p. and the value of 7. equal toy, say, 
for which 


Pop, ) = |», (an oetrl=a 
vives dlower | a contidence timit on the true (unknown) p. 


The limits p, and Pp, provide, of course, tworsided 1 2a confidence limits on the true 
(unknown), and in the tong run intervals so computed will capture p atleast 100C1 - 2a) percent 
ob the time, 


§ 


Thus, the tables of Harter? Leone, et al.? Hartley.” Thompson. and Vogler® may be used ta 


lind various confidence limits on the true unknown p of the binomial distribution. 


We remark also in passing that the percentage points for the binomial distribution are also 
confidence limits on the reliability of a “system” which consists of only one component for an 
observed number cof successes out of 77 trials in (esting the component or “system”. 


Example 2.14 
Was 30,¢= lO and a= .05, find solutions for Eqs 2-5 and 2-6. 
Solution 


Consider first the equation PU10,30.p) = .05. From the binomial probability tables herewith we 
read PC10.30.. 19) = 04508 and ?(10.30,.20) = 06109, from which linear interpolation yields p = 
1931. To use Harter’s table? we enter with 4 = ¢ = 10. Be a et 1 = 21, P= US, and read p = 
1933084. To use the table prepared by the Leone group.” we convert the equation to B(9,30p) = 
O85 and read p = 19330842. : 


For Eq. 2-6. we have B(1030.p) = .05 and in using the binomial probability tables it: is 
convenient to rewrite the equation as PU1,30.p) = 95. Then #C1130,.49) = 9384, POLL ,30,.50) = 
9500 and linear interpolation yields p = .4995. From the latter equation we observe that to use 
Harter’s table’ we need A = 11, B= 30 0 11 4 1 = 70, p = 98. But X€95.11,20) = 1 - 
VCOS20,11) = 1 - 5005613 = 4994387. From the table prepared by the Leone group? we read 
directly using the equation 1 its original ferm that p = .4994 3809, 
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It is interesting to note that since the Snedecor /distribudien Cand the Student &disGibutton) as 
related rather easily to the beta distribution, then binomial probability sums are expressible mi terns 
ot /-probabitities. By making a Cransformation in Eq. 2-1, it can be shown that 


| cent dna oan aanitaingean nr magmmaem me 


noe p : 
Boonpy = Pro] key ey) etree yoy 
etl lo op 

where /(r,. v2) follows Snedecor’s /-distribution with degrees of freedom. ep) = 2be + bband ys = 
Jar ¢). Setding Eq. 2-9 equal ta a yields . 
(ctl det 22a. 2c) 

ol Sa (10) } 
Qo ede ty gfle t Qan Qe) 


no 


° 





where the quantity 2) gl?) BP) is defined by. 

[hw db, oy 43) =| nN (2-14) 
The usefulness of Eq. 2:10 is somewhat restricted by insufficient (abulauion of / in terms ot 
particular degrees of freedom needed. However, the percentage points or standard probability levels 
for the ?"distribution are widely available. Nevertheless. it is not difficult to obtain binomial 
percentage points by previously discussed methods. 


3. THE NEGATIVE BINOMIAL DISTRIBUTION 


Consider a series of events of trigks for whieh Gi. Ce) and Ce) ef Conditions Cle} ) are satisfied 
but condition (da) is replaced by anew condition a’: 


(d') The series of experiments is performed a variable number of times until a fined number of 
successes (or defectives. ete.), say &, is achteved. 


Letting Y be the number of trials required to obtum exactly A successes, iC is easily shown That 
the probability that Y = y is 





| 
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As in the binomial probability case, we are usually interested ip a sum of probabilities such as Eq. 
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3-1. Denote the left-hand sum by : 
a Bt(kp) = > BX Ap) = PY Sed cak (2-2) 
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Then, we can show tat 

BY kp) = PkKepy (3-3) 
tn other words 

PY Se) =O X 2k) (3-4) 
where Y has a negative binomial distribution with constants k and p, and X has a binomial 


distribution with constants ¢ and p. To verify the truth of Eq. 3-3, assume that for a series of 
experiments, (a), Cb), and (c) of Conditions 1-1 are satisticd. Suppose we set as a goal the obtaining 
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of K successes but agree to perform no more than ¢ trials. “Then our goal can be achieved in two 
ways. Kirst, if the Ath success is obtained on either the Ath, A+ Ith. .... or the eth trial, the goal is 
accomplished. The probability that this series of events happens is PCY & o). ie., the left-hand side 
of Eq. 3-3. The second way to achieve the goal inc trials is to perform exactly c trials and obtain 4, 
A+ 1.2... c successes. The probability of the latter series of events is 2rCX¥ & &), the right side of 
fey. 353. dn both cases we are considering the event “the goal is achieved in c¢ trials” so that the (wo 
probabilities, therefore, must Le the same. 


If we define the right hand negative binomial sum by 


P*(c,k,p) 2 > b* (VAP) = PCY & e) (3-5) 


yee 
then we have immediately trom Eq. 3-3 


P*(aRK Py = 1 - PUKe ~ by) = B(K — le - bp) (3-6) 


Example 3-1 


Consider again the situation of Example 1-2. However, suppose now that sampling is continued 
until 5 detective items are observed, What is the probability that the Sth defective is found on or 


before the 40th item is sampled? 


Solution 


We have already discussed (a), (b), and (c) of Conditions 1-t, The new condition (d’) obviously 
applies here, Hence, we want Eq. 3-2 with & = 5, p = .08, and c = 40, We have from Eq, 3-3 and the 
tables 


B*(40,5,.08) = 2(5,40,.08) 


= .2132078 


Example 3-2 


To determine who pays for coffee, three people cach toss a coin and the odd man pays. If the 
coins all show heads or all show tails, they are tossed again. What is the probability that a decision is 
reached jn five repetitions or sooner? 


Solution 
A decision is reached on any given set of tosses if the result is two heads and a tail or two tails 
and a head. For any given set, the binomial is appropriate, and we identify 2 = 3 and p= .5 if a head 


is considered to be a success. Letting X be the number of heads appearing on the three coins, the 
probability of a decision is 


P(X = 1) + Pr(X = 2) 
= PAX > 1) P(X 3) = PC3,0.5) — P3,3.0.5) = 75 
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Concerning each set of tosses we observe 
(a) A decision is reached or not reached with each set. 


(bo) The probability of a decision is .75 for all sets. 


(c) The result of each set is independent of the result for any other set. 


(d’) A variable number of sets is required to produce a decision. 


in other words, the negative binomial conditions are satisfied with p = .75, k = 1. If Y is the 
number of trials required to produce one success (a decision), we want 


PrC¥ <5) = B*(5,1,.75) = PUL,S,.75) 


Using first Eq. [1-7 then Eq. 1-5 vields 
PULLS. 75) = B(4,5,.25) = | ~ P(S,5,.25) 
=} ~— .0009766 = .9990234 


If we fet Z denote the total number of failures obtained before the kth success is achieved, then Z 
+k = Y. Making this change in Eq. 3-1 yields the alternative form 


b*(z.kp)= (2 +4!) pk CL - py? 2 =0,1,2,... (3-7) 
k-1 


In some interesting applications (not arising in the context we have used), kK is not an integer. 
Then Eq. 3-7 is usually written as 
Me +k) 
z!IP(k) 


b*(z, kp) = pk (1) — p)?,z =0,1,2,. (3-8) 


Now, of course, Eq 3-3 no tonger appiies, and special tables are required to evaluate probabilities. 


~ Willian son and Bretherton” have prepared a table giving 
Ls 


ee Fla) = PZ <n) = >. b*(z,k.p) 
279 


to six decimal places for 4 wide range of p and k. The table seems to be quite adequate for a number 
of practical problems. For integer-vaiued k, the binomial tables are superior because of the fact that 


many more cnizies are available.’ 


4. BINOMIAL APPROXIMATION TO THE HYPERGEOMETRIC DISTRIBUTION 
In a number of sampling situations the following conditions are satisfied: 
(a) The drawing is made from a finite population of N items. 


(b) A random sample of size 7: is selected (drawing without replacement). 
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According to Eq. 4-5 
P(100,10,5,1) ~ BC,5,.10) = 1 — PQ,S,.10) 


Hence, 


piinetetie a ee ie re ee 


WXS 2) = P(2,5,.10) = .0814600 


Here, the actual value obtained from the Lieberman-Owen table is .076857. Of course, a higher 
degree of accuracy can be obtained with more complicated approximations. ' 


oe a 


hk 


Example 4-2 { 


a oe 


- yr ‘ 
pe 
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Approximate P(60,16,8,4). i 
Solution 


i According to Eq. 4-5, we have 


P(60,16,8,4) ~ a(48,38)- 1 -- P(S,8,.2667) 
Linear interpolation in the binomial probability table yields P(60,16,8,4) ~ .9641. In this case, the : 
correct value is 974228. : 


pe 


5. SINGLE SAMPLING BY ATTRIBUTES AND CORRESPONDING TESTS OF HYPOTHESES 


Suppose that a sample of » itens is selected at random from a lot and inspected, and that each 
item is classified as defective or nondefective (this is usually referred to as attribute sampling). 
Sampling is performed so as to make Conditions 1-1 realistic. The binomial distribution is hence i 
2 appropriate for computing probabilities for very large lots of items. Suppose that it is decided that : 
: if the true (unknown) fraction of defective items, say p, produced by the process is acceptably 
small, say Po (or tss), then satisfactory quality is being achieved, while for some p > py the quality 
is not considered satisfactory. Of course, in practice p will be unknown as indicated, and inferences 
will be based upon sampling. In order to make a decision we will observe X, the number of defective 
items in a random sample of size n, and conclude that the product is unsatisfactory if the observed 
value of X, say x, is too large. 





ee os 


Sometimes the sample size 2, the acceptable quality level py, and a small probability op 
(frequently .05 or .01) are specified. Then a rejection number d (so that c = d — } is called the 
acceptance number) is determined such that Pr(X 2 d) <p, if p = po. That is, we require 










Pld, n. pp) <A (5-1) 


In other words, if actually p = pp so that the quality is satisfactory, then we wish to reject the lot 
with low probability a). Since Eq. 5-1 will not yield a unique solution, we adopt the convention of 
making d as small as possivle. Observation in the binomial table immediately yields the required d. 
In terms of hypothesis testing we choose between the null and alternative hypotheses 


eS LAs ie OE NA Go 


Hy: P=Po and H,: p> po 
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selecting the critical region x 2 d so that the power is < QQ if p = po. 
Example 5-1 

Suppose n = 20, po = .05 and a,» = .05. Find the rejection number d. 


Sotution 





We seek the smallest d which satisfies P(d,20,.05) < .05. The binomial table yields P(3,29,.05) = 
.0754837 and P(4,20,.05) = .0159015. 
Hence, we sce that d = 4 (and c = 3). 


In sampling inspection practices and in many quality contro! type problems, it is frequently 
standard procedure to specify two values of p, say an acceptable level po and an unacceptabie level 
pP,, Where of course 95 < py, and to specify two small probabilities a and B,. If actually p = yo, 
then the product is regarded as satisfactory, and under such circumstances we would like to reach 
the correct conclusion with high probability, say at least 1 — ag. On the other hand, if p = p,, then 
the product quality is unsatisfactory and it is highly desirable that the probability of concluding 
incotrectly that p = pq be small, say no more than 8, . The problem is to determine the sample size 
n and rejection number d such that both conditions are satisfied. In symbolic form we need 


PAX 2 db <ag if p = po 
and 
PAX > d)>1 - B, ifp =p, 


In terms of &..0mial probability sums these inequalities are 


P(d,n, 7. 0 (5-2) 
Pd, > 21-8, (5-3) 
\ 
Usually, the minimum n is °-lred, at. sugh this need not be the case. In terms of hypothesis , 


testing we are choosing betwe 





Po and Ay: p=py ii 
by selecting a4 critical r-:gion x 2 a@ such that 
Power S a if p= po 
Power 2 1 — B, iip=p,; 

The small probability P(dnpy)) is the so-called “Type I error” or chance of rejecting the null 
ht thesis H_ when it is true, and 1 — P(d,np;) is the “Type II error” or the chance of accepting 
i hen actually the alternative hypothesis H, is true. 

If po and p, are entries in the binomial table, solutions for Eqs. 5-2 and ..3 are quickly found by 
inspeci.on. For a given d, the first inequality determines an integer 4, (d) such that the inequality is 


satisfied if n < n,(d). Similarly, the secone * quality is satisfied if m is large enough, say n > a 
n, (d). If n, (d) >n y{4) no solutions exist “- ‘e chosen d but if n, (d) <n y(4), then all » such ‘ 
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that n,(@@) <n<a yd) are solutions. It is easy to verify that as d is increased, both n, (d) and 
Ny(d) increase, but the upper bound a yp) increases at a faster rate. Hence, if solutions exist fora 
given value of d, they exist for all lurger d. If they do not exist for a given d, then they do not exist 
for any smaller d. Usually one would start with d = | and increase d a unit al a time until the 
desired solution is found. If after a couple of tries ii appears that d will be large. then d may be 
increased several units at a time until the vicinity of the desired solutions is reached. 


Example 5-2 

If po = .05, py =.20. & = .05, and B, = .10, find the sampling plan with minimum 7. 
Solution 

The inequalities, Eqs. 5-2 and 5-3 are then P(d,2,.05) < .05 and P(d.7,,20) > .90. From the table 
we find that with d = 1 we need 1 < 1 for the first inequality and n > 11 for the second. The 
essential sums from the table are 

P(1,1,.05) = .0500000 P(1,10,.20) = .8926258 

P(1,2,.05) = .0975000 PUAVY,.26) = .9141007 


so that with d = 1, we need 2 <1 and 7 = 11. Obviously, this is impossible so solutions do not exist 
for d = |. Similarly, we find that with 


2we needn < T7andn 218 
3 we needn < 16 and n 2 25 
4 we need n S 28 and 1 2 32 
5 we necdn < 40 and nn = 38 


RRR 


so that the desired plan is thus x = 38 and d= S (c = 4). Obviously, solutions for the inequalities 
exist for all d 2 5. For example, with d = 6, we find 45 <n < $3. 


Example 5-3 


If po = .10, p, = .20, ay = .05, and B, = .05, find the sampling plan with minimum 7. 


Solution 


The inequalities, Eqs. 5-2 and 5-3, are now P(d,n,.10) <.05 and P(d,n,.20) = .95. From the table 
we see that the first inequality cannot be satisfied with d = 1. Also with 


2weneedu< 3andza 222 
3 we needn < 8andn =30 
4 we needn < 14 and n 2 37 
5 we need 2 S 20 and n 2 44 


d= 
d= 
ad: 
d= 


Since it appears that it will take a while for the upper bound to reach rhe lower, we therefore omit 
some intermediate values of d. We find for 
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u 


d 10 we need i < SGandn > 76 
d iS we need 2h S 94 andn = 106 
d = 19 we need S127 anda = 129 
d = 20 we need n $135 anda 2 135 


u 


Hence, the desired solution is 2 = 135 and d = 20 (¢ = 19). 

When anv larger than the tabled entries is required and the p’s are not too large, we recommend 
using an approximation given by Hald.” Let x?7(v) have a Chi-square distribution with v degrees of 
. se oP 4 
freedom and define XV (v) by 

Pre OY < XO] = & 
Then Hald’s result yields 


4 l : 
wi y? (20 #2). aK : s)+e <n <— ]x? rer2). a ~ 5)te (5-4) 
9 1 vy Py 2 So Po 


To obtain solutions, ¢ is increased until the upper bound exceeds the lower bound and the interval 
contains at least one integer. 


Exampie 5-4 
Use Eq. 5-4 to solve Example 5-3. 
Selution 


From Faq. 3-4 we have 





3145 oe (2ce+2)+ (| <n <[95%0 (2e+2)+ «| 
By trial, we find 
with ¢ = 18, 49) 05-58) 418 Sng CaS 
or {29.11 <n < 127.18 
and with ¢= 19, G5)S5.79) 119. <y < £9.5)(26.51) +19 
or 134.96<n < 135.42 


Thus, the approximate procedure also yields 
n= 135 and c= 19, Chisquare values were obtained from Harter's? table. 


The method illustrated in Exampies 5-2 and §-3 is given in a paper by William C. Guenther, “Use 
of the Binomial. Hypergeometric and Poisson Tables to Obtain Sampling Plans”, Journal of Quality 
Technology, Vol. t, pp. 105-109, 1969. 


XVI 
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In an cartier reference)? Grubbs (1949) gave a solution to the problem of oe the single 
sampling inspection plan. a and oc, such that as operating characteristic curve (i. e.. one minus 
power) passes very near (wo points (1 yo) and (8, .p,) for risks of a) = .05 | B, = .10. In 
this connection. one has only to find frem the ay = .05 table those sampling inspection plans or 
various values of, and ¢, along with o's very near the “desired” po for which P(c,21;po) = 95, 
and corresponding plins 7, and c, with p,’s near the “undesired™ py for which P(c,aejp1) & .10. 
Then, out of the two sets of plans. there is a common one, 2 and c, which guarantees the desired 
risks (very nearly) and for which a is the smallest sample size meeting with the stated conditions. 


DISTRIBUTION-FREE 6-CONTENT TOLERANCE INTERVALS 


Let Vy. Vy..... 4, be a random sample from a continuous distribution having unknown density 
function /(v) > O. where a <x <b, and a and b are known. Denote the ordered sample by X( 1). 
Nosy... X,, and let 


(n) 


an 

= PAX, <X SX) =f 1 fy) de i<j 
x 
(0 


Wis thus the fraction of the distribution captured between the ith and jih order a Denoting 
the observed values of AZ, and Xun by x,,, and X,,,. respectively, we call (x,y. ¥ a two-sided 
tolerance interval for the distribution of Y. The numbers Xj) and x yy are usually Baa to as the 
lower and. upper tolerance limits. respectively. The intervals (a, x,,.)) and (x b) are called 


; : ) 
one-sided tolerance intervals. 


(1)" 


Suppose now that we are given B and y, two positive numbers less than i (but usually close to 1). 
Wand / are determined so that 


Mii 
lr = B) = olf fix) dv 2B |2y¥ (6-1) 
, 


omar) 


then (v,;)- Xyy is called a B-content tolerance interval. It can be shown that Eq. 6-1 ean be 


expressed as 
Pr ftitidoni -p2y. (6-2) 


a formula which also holds for the one-sided cases if we agree to set ¢ = O when seeking only an 
upper tok tance limit and to set / =a + 1 when seeking only a lower tolerance limit. Since Pr(W 2 B) 
does not depend upon f(x), the tolerance intervals are “distribution free” 


Sometimes (as with sampling plans and hypothesis testing) four constants By <8,, and 79 > 7, 
are given, It is then desired to obtain f and / so that 


Pr(W = Bo) = Yo (6-3) 
and 
MrW2B,)<y, (6-4) 


In descriptive terms, it is desired that the tolerance interval should be long enough to capture (with 
high probability) at least a given proportion of the distribution, but also that it should not be 
unnecessarily long. The second condition is met by requiring the interval to capture a iarge 
proportion B, of the distribution with low probability 7,. We can replace Eqs. 6-3 and 6-4 by 


xvii 








Wade abd Bade el ek 


a 







SRI aE ra emecmeR meee 


ae or : 


AMCP 706-109 


Pn —j tit lot —B d= Yo (6-5) 
and 
Pin -jritdad—-B)<Syn (6-6) 


The three main problems considered with distribution-free B-content tolerance intervals are: (1) 
Given n, B, and y, find i and j (one-sided cases), or } — i (two-sided case) so that Eq. 6-2 is satisfied, 
and Pr(W = B) is as close to 7 as possible. (2) “iven /, 7, B, and , find the minimum ” which satisfies 
Eq. 6-2. (3) Given Bo, By, Yo, and y;, fin. av ininimum ¢# and accompanying & j, or j — ¢ which 
satisfies both Eq. 6-3 and Eq. 6-4. Although various tabulations fer these problems have appeared in 
the literature, all three can be solved by inspection with the probability table. 


oT eter Shere a ed 





S Example 6-1 


If 6B = .90, y = .95, and n = 50, find the distribution-free one-sided and two-sided 6-content 


tolerance intervals with Pr(lW > .90) > .95 being as nearly satisfied as possible. [f 6 = .90, 7 = .95,7= 
!, and j = » — 1, find the distribution-free one-sided and two-sided £-content tolerance intervals 





with minimum sample size. 





Solution 
For the first problem, Eq. 6-2 yields for each of the three cases 


P(i,S50,.10) > 95, P(S1 — j,50,.10) > .95, and P(S1 — j +7,50,.10) > .95, 


By inspection we find i = 2, Sf — j = 2 or j = 49, and $1 -- G — i) =2 orj — i = 49. Then, for the 
one-sided cases the desired intervals are (cy ~) and (=, Xag y)- For the two sided case we have 


Way soy) 
For the second problem, Eq. 6-2 gives 


P(,n,.10) > 9S, P(2,7,.10) = 95 and P(3,n,.10) > .95. 


ow Lat al $i Fata ae ellie ad rt yiseleer iy Z, iabatal hea Nes ei pin a or ~re b Poe a 


By observation in the table we find that to satisfy the inequalities we must have n > 29, n > 46, and 
n > Ol, respectively. Then, the desired intervals are (x,,,, ©) obtained with w= 29, (—», x ) 

; 2 : ; aA , (45) 
obtained with n = 46, and (x, 1)? % 6 0)) obtained with a — 61. 


or | a 


Example 6-2 
If By = .90, B, = .99, yo = .95, and yy = .05, find the distribution-free one-sided and two-sided 
Bo-content tolerance interval with minimum 2. 


had a 


me Sts 


Hehe S Solution 


Consider the two-sided case. Eqs. 6-5 and 6-6 are then ne 
. 13 
‘ Pn — j +i + 1,n,.10) 2 95 7 
“3 
: and Pin —j+it 1.01) <.05 3 
x 

The solution is obtained by trial, starting with n — 7 +i+ 1 = 1 and increasing this quantity until a Es 


solution is found. With 1 — j++ 1= 1, we find by observation from the table that we must have n 
> 29 and n <5 to satisfy the two inequalities so that no solution is possible. With n —- jf +it 1 = 2, 


. 
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we need a = 46 and a & 35, again an impossibility, Withn fj +it+ f= 3 we needa Sol ands 
82. The desired is thus 61 andj - i = 59. The tolerance interval can be either (xo, 6X6 4 )) OF 
iy Xeon) 


With a lower tolerance limit we need /= 3 and 2 = 61, the tolerance interval being (x, ,,. 7). With 
an upper tolerance limit, we need a -j +7) +1 3. 2 61, and/ = 59, so that the interval is (+. 


X¢so)) 


The material of this paragraph is taken from “Determination of Sample Size for Distribution-Free 
Tolerance Limits”, William C. Guenther, The clinerican Statistician, Vol. 24, 1970, pp. 44-46. 


7. DOUBLE SAMPLING PLANS WITH TWO COND! TIONS ON POWER (OR OC)* 
In par. 5S, we considered single sampling plans which satisfied the two conditions 


Power Sa if p= pry 
(7-1) 
Power 21 > B, ip=p, >» 


and observed that in the equivalent hypothesis testing problem we chocse between 
Ho: PFPo and Wyo pep, 


Now we discuss the determination of double sampling plains which satisfy the same iwo 
requirements. The type of double sampling scheme which we consider here requires two sample 
sizes a, and Ay. and two rejection numbers a, and d@; (or two acceptance numbers ¢, Fy -~ Lard 
Cy =a, ~ 1). The procedure begins by inspecting a sample of size ny items and counting the 
number of defectives (or successes), say VY). Ifxy Sc,, where x, is the observed value of X,. we 
accept On the basis of the one initial sample. Simiiarly, ifs, > c2, we reject without (aking a second 
sample. However, if ¢, <x, <cy, a second sample of size ay is taken and the nuinber of defectives 
X, observed in the second sample recorded. Then. we accept ify, +42 Sey and reject ify tx. > 
¢). The major advantage of doubie sampling is that on the average fewer observations are required 
to reach a decision. (In “screening” or “rectifying” inspection - sce par. 9 on Dodge-Romig LTPD 
Plans - less average total inspection is necessary.) 


There will be many sampling plans 24, 7y. ¢;. and cy satisfying the given conditions. (in fact. as 
in the single sample case. there are an infinite number of solutions.) Consequently, in order to 
obtain a unique plan we can impose an additional condition. Actuaily, this was also done in the 
single sample case of par. 5, where we minimized the function GQz.c) = 2 subject to the given power 
requirements. In double sampling we need to consider a function of four constants 2). 2, ¢,. and 
Gy, say GCM, Mz. Cy. C2). Possible candidates for such a function inchide 


(a) Any function of the average number of observations required to reach a decision. 


(b) The Dodge-Remig!! (and Eq. 9-2 that follows) average total inspection number required in 
screening or rectifying inspection. 


(c) A generalization of the Dodge-Romig function to a cost function as proposed by A. Hald.!? 





*Operating Characteristic. 
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It is the responsibility of the individual designing the sampling plan to construct his own G-function 
so that it reflects some criteria of interest. The function can be minimized (or maximized) by trial. 
In this paragraph, we consider only plans which satisfy Conditions 7-1 and not the minimization (or 
maximization) of a G-function. 


For the type of double sampling plan which we have described, it can be verified that 


Power = K(p,m,,2,€,.d2) 


d,-d,-\ (7-2) 
= P(dz.nyp) + S bd, +iny.p) Pd, — dy —iin2p) 
' 720 
Consequently, we may write the two conditions of Eq. 7-1 as 
K(pu.my hy ay d+) < Xo (7-3) 
and 
K(py itty ta. Gy ,d2) > l = B, (7-4) 


Obtaining solutions for Eqs. 7-3 and 7-4 is accomplished by trial. If we use the information that 
follows, these solutions can be found rather rapidly. 


To determine which double sampling plans satisfy the given power conditions, we begin by listing 
the single sample solutions (and nonsolutions). This information is essintial to the double sample 
case. As demonstrated in par. §, we can obtain the following type of information: If 


d=) nSNy(Q), 2 en, (1), 2, QD) > ny, no solution possible 


d=2 nSNy(2),n Sn, (2), a, (2) > ny QQ), no solution possible 


d=r--} nSny(r— ln Fn (r- 1), nr D> ny ~ 1), no solution possible 
d=r n<ny(r),n Zn, (0), solutions for ry) <n <ny) 
d=r+l nSnyti)nen, et 1), solutions for n@ttsne nyt {) ete. 


That is, solutions exist if d >r but not if d<r. Moreover, for the double sample case, the following 
comments are helpful: 


(1) Double sampling plans satisfying the given power conditions exist only for alld, > rand d, 
= 0, 1, 2,..., d, — 1. If solutions exist for given d, and @, with n, = 1... the minimum r to 
admit solutions, then solutions exist for all 2, such that #4, <1, </,(d). a set which contains 
the single sample solution. Ordinarily one would be interested only in solutions obtained with d = 


r, since larger d's require larger sample sizes. 


(2) Bounds on n, for chosen d, and @, are immediately obtainable. It can be shown that we 
must have 


Pd, yp) e 1 By (7-5) 
which means ny >myz, and 
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P(dy.8y. Po) SX (7-6) 
which means ny Stru- 


Solutions do not always exist for a, = yz. but this number provides a lower bound in our 
search. 


(3) It can be shown that a, +”, 2 2, (d2). In other words, the combined sample for double 
sampling has to be at least as large as the minimum sample size possible in single sampling, 

(4) For given d, and d,, lower and upper bounds on 72, say az, and #,y, can be found for 
cach n,. Of course, if mn. < 227, solutions are not possible for the chosen ,, and #, has to be 
increased. Both 23, and my are nonincreasing functions of n, which behave very well, a fact 


which considerably reduces the number of trials required, since as ay increases 2, + 27 tapidly 
approaches n, (d2) and equals 1, (d>) for at least ny =", (d2). 


In view of these observations, the recommended procedure is 
(1) List the single sample solutions and nonsolutions. 
(2) Select any d, (usually d. =r) for which solutions exist. 


(3) Select any d, «ach that | <d, <d,. We would probably begin with d, = 1 and increase d, a 
unit at a time. In a number of plans used in practical situations we have d, < 4, /2. 


(4) For chosen d, and d. determine bounds on »,. We would probably begin looking for 
solutions with the smallest 2, satisfying Eq. 7-5 and increase 2, a unit at a time. 


(5) By trial for the chosen @,, dz and n,, find bounds on n, determined so that the two 
conditions on power are satisfied. If 22, <n2y, then one set of solutions is the chosen @,, d,, 4. 
and nap <n, Suzy. Usually, we would prefer that 2 be as small as possible, ie... = N27. 


(6) Select another n,, d,, @; and repeat Step (5), proceeding in some orderly fashion. In an 
application we would terminate the search when we have found the solution minimizing (or 


maximizing) the appropriate G(r), ny. C1, ¢2). 
We now consider an example. 
Example 7-1 
Find some double sampling plans for the case y»9 = .05, py =.20, & = .05, and B, = 10. 


Solution 


In Example 5-2 we considered the sing!e sample problem using the inequalities (d,n,.05)< .05 
and P(d,n,.29) > .90 


d=1 n& land ntl, no solution possible 
d=2 n<& Tand n> 18, no solution possible 
d=3 nS l6and nn? 25, no solution possible 
d=4 


with 


nS 28 and no 32, no solution possible 


; 
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d= ws i ness, or 38 <n < 40 
v=o ” neds. ordsS <a S S83 


ele, 


Double stmpling plans are thon possible for dy # 5. With dy = 5 we can have dy = tL, 203. or 4. 
Since all solutions are oblaincd ina Similar manner, we consider only the cuse dy = 4 and, = 5, 


lneguibties, Eqs. 7-3 and 76, are now P4yny. 620) 2 .90 and MS, .05) < 05. From the 
previous single sample results we see dat these (wo inequalities imply 2, 2 32 and n, < 40. 


Withe, = 32), =4 and dd, = 5. the power at p = .20 is 
Kb. 20.32.01 4.5) = POS 32. 20) + 0(4.32,.20) PU ty 20) 
= 79562 + 11129 Plan. 20) 


Using the fact that ay ter. > UR we find by trial 


K 


.20,.32,13.4.8) = 79562 + ET E29¢.94502) = .90079 
and KC.20,32,1.2.4,.5) = 79562 + 11 129¢.93 128) = 89926 


so Ura C2 makes power 2 OG. Similarly, with p = .05 we have 


K(OS.32 203 4.5) = PUS,32, 05) + b(4,32,,05) PUL ar 05). 
Ky trial we ootain 
K¢.05,32,15,4,5) = .04904 and A(.05.32,16.4,5) = 05028 
so that W. < 15 makes power & .05. 


daving found solutions with one 2,, less labor is required for a, + 1. With a, = 33 se 


KO.20,33,8.4.5) = 90219 and AC 20,33,7.4.5) = 89794 so that uz. 2 8. Also. K(.05.33.12 
04959 Kt.05,33,13,4,5) = O5115 so that nr, S12. 


m4 
vs — 


Contig the calculations. we soon have all the solutions ford, = 4 and d. = §. These are 


ry? 3d, 13 <n, <5 = 37, }<n, <3 
ny 733, 8<Sn, S12 i, 58. O<n, <2 
n, = 34, 6S, 9 ny = 39, O0<n, <1 
ny, = 35, 4<n,S 7 ny = 40, n, = 0 
= 36 3<n, S 5 


Similacly. we could obtala all the solutions for d, 1, 2, 3 with d, = §. Had we sought to 
minimize a G-function, it would have been better to start with d, = 1 followed by d, = 2, since d, 
= 4 yoquires such a large ay 


Since the secon¢ sample may or may not be taken, it is of interest to know the average number 





4, 
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J 
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‘ at observations which are required to reach a decision (ie.. tie ASN of average sample number) 
with any double sampling plan a2). ay. d,. and dy. We see thats, observations are required i vy S 
7 ay Lor vy 2 dy. while a, + ay observations Must be taken fd) ts vy <a a. Hence, the random 3 
i vartable Iv, the total sampic sive. takes on the value a, with probability Pray, <d,o- 1) + PAN, > + 
an - d,) and the value a, toy with probability Pra, < VY, < a2), and by the definition of expected = 
| value We have that? FE 
od EUW) = ny LPP dy Sy — DE PN, Bay) ty tag Prd SX, <d2)| ¢ 
N ig u 
. zany tay lPrd, SN, Say) q 
: { a 
In terms of binomial sums we have then that ECW) is given simply by 3 


ECW) sony tary (Pld typ) Pda. 2) (7-7) 


ae 
Oren ey 


tor any quality level p. 


om 


ld 


8 CURTAILED AND SEMI-CURTAILED ATTRIBUTE SAMPLING: SINGLE AND DOUBLE 
SAMPLE CASES 






Suppose that a single sample attribute plan is vo be used. Sampiing is performed so that the 
Binomial Conditions 1-1 are satisfied. The sample size mand the acceptance number c¢ (or rejection 
number ¢ + 1 =d) have been determined possibly to sausty two conditions on power as in Example 
5-2. The usual way to proceed is to take u observations. observe the number of defectives VY. and 
accept if and only ifx <¢. We will call this procedure regular binomial sampling. 





lth _ eer, 


Eee 


If observations are expensive. there are a number of ways to attempt to reduce the amount of 
inspection. One of these, called curtailed sampling, calls for taking observations one at a time and 
quitting as soon as the decision is obvious. To illustrate, consider the sampling plan of Example 5-2 
4 where we found a = 38 and c = 4, If items are inspected one ata time, it is obvious that we could 
stop sampling and conclude that the process is unsatisfactory as soon as the Sth defective is 
observed. This could happen as early as the Sth observation or as late as the 38th observation. 
Similarly, once we have observed 34 nondetectives, then sampling can be terminated with the 
conclision that the process is satisfactory since it is then impossible to obtain more than 4 
defectives in 38 observations. Since ¢ is usually rather small compared to a, a reasonable 
approximation to this curtailed sampling procedure is obtained by stopping as soon as too many 
defectives are observed. When this is done, the term semi-curtailed sampling is Used. Intuitively, we 
would expect the power or OC curve to be the same whether we use regular binomial sampling, 
curtailed binomial sampling, or semi-curtailed binomial since the latter two procedures can be 
regarded as special cases of the first in which the inspector does not record all the information 
available to him. That our intuition is correct can be verified mathematically. In general, with 
curtailed sampling we stop sampling and (a) reject as soon as c + 1 defectives are observed, or (b) 
accept as soonas a - ¢ nondefectives are observed. In the semi-curtailed case We reject as soon as ¢ 
+ 1 defectives are observed, accepting if this fails to happen (in which case all observations are 
made), 


ee Oe 


| 
| 
: 
| 
| 


When sampling is curtailed or semi-curtailed, the number of observations actually required to 
reach a decision, say W, is a random variable which can assume values frome + | to a. It is of 
interest to Know the expected value of W, say ECW), sometimes called the average sample number 
(abbreviated ASN). Obviously. this average will be Jess than # and will depend on p. If sampling is 
curtailed, it can be shown that 
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‘ ; c+] noe : 

{ Law) = se Pet Sonet i PY] ttl Met haat typ] (8-1) 
* 

; The corresponding formula for scini-curtailed sampling is 


sh 
L(y) = Me tQat lp) bath - Pet dead (8:2) 


nlm 


%Z 


Example 8-1 


For the sampling plain of Example 5-2 find the average number of observations required for yp = 
OS. 2O..18 and 20 with both curtailed and semi-curtailed sampling. 


“7 


ons tp BD ae al! ate 


Solution 


We hada = 38 and ¢ = 4. For curtailed sampling we need 


ee OG Ad thes 
EUW) serra LASCSRIAL acer (l PCS.39p)) 


With p = 05, we have 


EW) = “Fe (.01 230) ge A370) = 38.45 


wat ad eS, ata eo 


eae 


Similarly, with p = 10. 15. and .70. we find 


\ 
by 





LOW) = 34,05, 29.43, and 24.18, respectively 


ee Dl ate ae 


4 When sampling is semi-curtailed, we need 


Pia ne 


KW) = 416.9991 $38 (1 - PU38p)1 


With p= 08. we obtain 


ECW) = Zz (.01230) + 38 [1.03973] = 37.72 


Similany, with p= 10. 15. and .20, we tind 


EW) = 34.98, 29.70, and 24.25, respectively 


Curtailed and semi-curtailed sampling can be used in conjunction with double samrpting to further 
reduce the expected number of observations. Formulas tor the ASN. however, are niuch more 
complicated than they are for regular double sampling. Since there is hitthe difference in toe ASN's 
obtained from curtailed and semi-curtailed sampling. wd formulas for the latter are slnpler, we will 
consider only the case in which scmi-curtailed samp! ag is used in both samples. In Chis case then, it 
can be show that 


eee ee er ee 


eevee” 


Of 
eae Mm ae ae 








AMCP 706. 109 


er ds 
haw) ape Pldyt Lang t Lad tary | Pldy ay py 


4 
: t Hy Pd, vty y) Kipary ty dy a; | se) 
fe ahy 11 
biog S (dy dy OC HA Ptdy dy HL fats BY) 
D 
bow 
270 


where AQpar vty dy ds) is defined in pars 7 by bq. 7-2. As in dhe single sample case. the power is 
the same for regular, curtailed. and semi-curtailed double sampling. Cran! * gives a result equivalent 
to Eq. 8-3 which can be used when poy. and ary are not entries in the table, 





Tae ee ee 


| 
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Example 8-2 


seer 


If semi-curtailed sampling ts uscd with both samples tor the double sampling plana, ~ SQ, = 
1O0,d, # and dy = 4 find the ASN for p = 05 and p *.10. 


4 scananil : 


Solution 


4 


aH \ 
“4. 


f 
td oust flab 


We use Eq. 8°3. With p = 05, the tirst Qo terins are 


ogy ult, ya 


i tae 


4 oe ee biuss 
OE POSTS) t+ SOql  P(4,50..05)) 


ay 





LUN A eM, 


= SOC TOA) | S00 Tod) = 40.8537 


cv avid a it oa 


For the next term we peed 


PEVSO.OS) + bE2,80, 05) POL 10U,.08) 63,580,058) POTO0,.08) 


it 


K(.05,50,100.2,4) 


tee + 
1 


t 


= 239892 + C2610) C9GQYT9) + TPIS TS) COOdO8OD = .7OOSS4 
Then the thitd term is 


: 1001(2,30..05)  K(.08,80.100.2.4)] 
| 
= {Q0].720568 700584] = 1.0984 | 
The last ternr is 
oe hgl 20250 08) PQ3, 101.05) -# b(3.S0,.05) P(2,101,.05)] = 


Finally, we have EUW) = 46.8537 + 1.0984 + 13.4926 = 61.44 


Similarly, when p = .10, £(W/) = 35.93. 


By Eg, 7-7. it can be verified that with regular double sampling the two ASN’'s are, respectively 
98,10 and 71.65, which are therefore much larger. 


i 
20{ 20.261 101) C.885796) + (219875) (.964477)] = 13.4926 : 
' 
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9. DODGE-ROMIG LTPD (LOT TOLERANCE PERCENT DEFECTIVE) PLANS 


Consider a product that is mass produced and is assembled at random into lots of size NV. From 
each Jot, items are drawn at random and the number of defectives in a sample of n items is 
observed. If the lot is accepted on the basis of the sample of 7 items, then all defective items found 
when sampling afe replaced by nondefectives. If the lot is rejected on the basis of the sample, then 
the remainder of the lot (N -- 1 items) is also inspected and all defective items found dering 
inspection of the lot are replaced. The procedure just described is called screening* or rectit'ying** 
inspection. Our goal here is to find single and double sample rectifying inspection plans which 
minimize the average total amount of inspection if the process fraction defective p is equal to 77, the 
“process average”, buf subject also to the condition that the powerbe> | — 6, (or OC < (, ) if the 
process actually produces a fraction defective cqual to an undesirable level p,. (in the language of 
Dodge-Romig", p, = Pp, = LPPD. Le., the tot tolerance percent detective, and B, = the Consumer's 
Risk.) 


We consider first the single sample case. A sample of size n is selected at random and the number 
of defectives Y is observed. If x < d defectives are found in the sample, then the lot is accepted 
(defective items in the sample replaced), while if . 2 d the lot is rejected and all defective items in 
the lot replaced. The unconditional distribution of X is a binomial probability distribution with 
paraineters 7 and p, and the power is clearly given vy Pr(X > d) = P(d, np). The average tmount of 
inspection per lot is thus seen to be 


ig Fat (N ~ n)P(d,n,p) (9.1) 


An LUPD rectifying inspection is obtained by minimizing Eq. 9-1 when p = p subject to the 
condition 


P(dvayp,) 2b -- By (9-219 
(To be strictly correct, the left-hand side of Eq. 9-2 should be a sum of hypersvometric 
probabilities. However, as Dodge and Romig'! have done. we use a binomial approximation.) The 
minimization is accomplished by trial by starting with d = | and increasing @ one unit at a time. For 
each d, the minimum 7 satisfying iq. 9-2 is found by inspection in the table, and /, iscompuir d for 
p =p. Calculations should cease of course when the minimum, say I,, is observed. 
Nodge-Romig!! give some plans of this type for 8, = .10 and for a few values of p,. 


Example 9-1 


If iv = 100, py = .10, and 8, = .10, find the rectifying inspection plan which minimizes the 
average total inspection for an average process Jevel of p = .02. 


Solution 
According to Eq. 9-2, we need P(d. n, 10} 2 .80. Using the binomial table, we find that with 
d = Landn = 22,1, = 22 + 78 P(1,22,.02) = 22 + 78(.35883) = 49.99 


7A term due to Dodge and Romig. 
**A term due to the British. 


XXVI 
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iy 
Kt 


2 and n = 38, I, = 38 + 62 P(2,38,.02) = 38 + 62(.17603) = 48.91 


d 


Mi 


3 andn = 52, I, = 52 +48 P(3,52,.02) = 52 + 48(,.08593) = 56.1) 


further calculations are obviously unnecessary. The desired sampling plan is therefore nu = 38 and d 
= 2 (c= 1) with/, = 48.91. 


Next. We consider the somewhat more complicated double sample case. From each lot, a random 
sample of sizer, is drawn and inspected. Let X, be the number of defectives in the sample of a,. Uf 
x, <d, defectives are found in this sample. the defective items in the sample are replaced by 
nendelectives and the fot is accepted without further inspection. ff.<, “2 d,, the Jot is rejected 
and totally inspected, and all defective items are replaced by good ones. If dé, Sy, <2. a second 
sample of size vy is selected at random from the remaining N~ ay itens. and ¥,, the number of 
defectives in the second sample is observed. Then. if, Sx, tx. < dy. the fot is accepted without 
further inspection, bul all defective items in both samples are replaced by good ones. Finally. if x. 
+ x, 2d, then the lot is rejected and all defective items in the lot are replaced by acceptable iterns. 
Tire unconditional distributions of Y, and X, are binomial with parameters p and ,, and p and ny. 
respectively. Now the average total inspection for the double sampling case is 


Spel bie at at eet ee ae 


Ry. 


Py ty ty PCdy ay py t(N ony tg) K(pyiny ey .dy dz} (9-3) 


Mens Maske neces Ree. 


where the K-function is the power defined by Eq. 7-2. We seek to minimize Eq. 9-3 when p= p, 
getting say /),, subject to the conditon 


ot 


jae 


K(py ont) ats .dy.dy) 21 - By (9-4) 


Aga. minimization may be accomplished by trial. (As in the single sample case, to be strictly 
correct the left-hand side of Eq. 9-4 should be given in terms of hypergeometric probabilities.) 


wm erate 


pagers 
NX 


Since there are now so many possibilities. we formulate a systematic approach. 
The following comments are helpful: 


(1) For given a,,d,, and d,, select #3 as small as possible so as to satisfy Eq. 9-4. A Jareet a, 
willjust make /,, larger. 


> 
(2) Kor every d, end d;, the condition OC < 6, for p =p, implics 
PMdy ny Pye b- By (9-5) 
This inequality yields a lower bound on 1, , giving 7, 2 4 ),. say. Also if 
Pld, ty .pry2) — By (9-6) 
then the power >| B, for every ny. Letting 2, be the smallest 2, to satisfy bq. 9-6, we need 


consider only #, <a, co. Usually it will be apparent that calculations with the given d, and d2 can 
cease before n, is made as large ag yy. 


(3; Only for 7, < /,do we need consider the plan. since otherwise 2,, 2 fo ane the ebieetive of 
dcouble sampling, to reduce average total inspection, would not be possibic. 


& 
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(4) For chosen. d, and d, we must haven, tn, 2 ayy. 


(S) For chosen d, and d@,, it may be shown that as », increases, the sum n, + ny is 
nonincreasing and has a minimum yalue of 7, y. 





In a numerical problem, we suggest the follow’. steps: 





7 4 (1) Calculate 7, since, as we have already mentioned, it is not necessary to consider plans for £ 
a which /, > [,. ) 4 
: a (2) With d, = 1, determine ,; from Eq. 9-5. 4 
z ay ' - == 
: (a) With dz =2 F 
ss . i F 
¢ 1. Find n,y from Eq. 9-6. q 

: Fl 


2. By trial, using the fact that a, + 2, 2 ”,y, then find the mmimum value of ”,, say : 


Ny 1, Which satisfies Eq. 9-4. 


sal AI bans et 


3. Withn,, andn2,, find Z, for the process average p = p. ; 


th 


(b) Repeat (a) with d, = 3. As a first guess, for the new mz, increase the old na, by the 
same amount , y has increased. 


(c) Repeat (a) with d, = 4, etc. 
Terminate the process when it 1s cbvious that /, must increase with any further increase in 
dy. 


(3) Repeat Step (2) with a, repiaced by n,,, + 1. Then repeat Step (2) with 2,1, replaced by n,, 
+ 2, etc., terminating the process when it is obvious that a minimum has been found for each d, and 


for which it has heen necessary to consider d; = 1. 


(4) Repeat Steps (2) and (3) now with d, = 2, then with d, = 3, then with d, = 4, etc., and 
terminate these trials when the values of J, become too large. This happens at worst when #1, > Is, 


maybe sooner. 


(5) By observation, select Z,. 


The major advantage of these suggested proceedings is that in the calculation of power all the 
previous b(d, +/,1,,p) are used asd_ is increased by a unit. We now consider an example. 


Example 9-2 
If N = 100, py =.10, and B, = .10, find the double sampling rectifying inspection plan which 
minimizes the average total inspection for p = .02. 


t 
[ 
t 
Solution | 
In Example 9-1, we found T,= 48.91. Hence, we need consider only those plans for which J, < | 


48.91. 


XXviil 
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As a preliminary calculation, we find the minimum n to satisfy P(d.n,.10) > .90 for a few d. 
From the binomial table, we find for d = 1 that n = 22, ford = 2 that n = 38, ford = 3 that » = 52, 
for d = 4 that v = 65, ford = 5 that 2 = 78, etc. Because of Eq. 9-5, the only candidates ford, are 
seen to bed, = 1 and d, = 2 since with d, > 3,n, 252 and dp) > S2. 


With d = 1, we need according to Eqs. 9-5 and 9-6 that 22 <n, < 38. Also, we must have 7, + 
n, > 38. With d, = 1, d, = 2 and, = 22, the power at p = .10 is 


K(.10,22,12 1,2) = P(2,22,.10) + 6(1,22,.10) PCL,n2 ,.10) 
= 66080 + .24072 P( 1,3 ,.10) 
By trial, we find K(.10,22,48, 1,2) = .89999 and K(.10,22,49,1,2) = .90014. 
Then with 2, = 22,2. =49,d, =1, and d, = 2, we have 
K(.02,22,49,1,2) = P(2,22,.02) + b(1.22,.02) PU ,49,.02) 
= 07096 + (.28787) (.62840) = .25186 


Thus, 


I = 22 + 49 PUL,22,.02) + 29 K(.02,22,49,1,2) 
= 22 + 49(.35883) + 29(.25186) = 46.89 
We next take d, = 3 withd, = 1 anda, = 22. 


Now n, tot, 2 52, since 7, y = 52. When p = .10, the power is 


K(.10,22,14 1,3) = P(3,22,.10) + b(1,22,.10) P(2 722 ,.10) 
+ b(2,22,.10) PC, ,.10) 


= 37996 + .24072 P(2,n3,.10) 
+ 28084 P(1,n,.10) 


By trial, we find 
K(.10,22,69,1,3) = .89987 and K(.10,22,70,1,3) = .90002 
With #2, = 22, n, = 70,d, = 1, and d, = 3, we have 


K(.02,22,70, 1,3) = P(3,22,.02) + (1,22,.02) P(2,70,.02) 
+ b(2,22,.02) P(1,70,.02) 


= 00927 + .28787(.40956) + .06169(.75688) = .17386 
and when p = .02, we find that 


Ty = 22 + 70 P(L,22,.02) + 8K(.02,22,70,1,3) 
= 22 + 7.35883) + 8(.17386) = 48.51 
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Since the second term of /, is increasing faster than the third term is decreasing, it is not necessary 
therefore to try d,; = | and d, = 4. 


We next cepeat the above calculations with n, = 23, then with n, = 24, then with 2, = 25, and 
ivrminate the calculations when it is obvious that minimum /, have been found for d, = t, d, = 2, 
and for d, = 1 and d, = 3. The results are as follows: 





d,=t,d,=2 d, =1,d,=3 
(my.n2) Ty (n1,m2) Ip 
(22,49) 46.89 (22,70) 48.51 
(23,29) 43.74 (23,48) 44.47 
(24, 23) 43.16 (24,41) 43.65 
(25, 19) 43.00 (25,37) 43.66 
(26, 16) 43.06 (26,33) 43.50 
(27,14) 43,49 (27,31) 44.08 
(28,12) 43.82 (28, 29) 44.60 


Fiaatly, we repeat all the above calculations with d,; = 2. Scame of the results are 
i 





d, =2,d, =3 d,=2,d,=4 
(ny, tty) Ty (t,,72) Ip 
(38,34) 47.00 (38, 54) 48.06 
(39, 24) 46.98 (39, 53) 47.92 
(40, 20) 47.56 (40,38) 48.42 
{41, 16) 48.03 (41,33) 48.81 
(42, 14) 48.80 (42, 30) 49.48 


We observe finally that the desired plan is ny = 25, n, = 19,d, = Land d, = 2 with J, = 43.00. 


Since B,; = .10 and p, = .10, we can compare our selection with the results tabulated by Dodge 
and Romig! !. They give # = 20 and d = | for the single sampling inspection plan. Since P(1,20,.10) 

= .87842, the power requirement is not quite met. This plan gives J, = 20 + 80P(1.20,.02) = 20 + 
80(. 33239) = 46.59. Recall, however, that we had # = 38, d =2 and I= 48.91. For the double 
sample case, Dodge and Romig give mn, = 25, n, = 13. d, = 1 and dz = *y The power of this plan 
when p = 10 is K(.10,25,13,1,2) = .87752, again not quite meeting the specified requirement. The 
plan also yields J, = 40.06 with p= .02. The Dodge-Romig solutions were tabulated for intervals on 
both N and p. and in most cases the ”’s were rounded to the nearest 0 or 5. Hence, both tabled 
plans are reasonable approximations to the plans we found. The main advantages of the procedure 
{ of this paragraph over the Dodge-Romig tables are: (a) plans exactly meeting the specified 
‘ conditions can be found, and (b) 6, can be any value or level we choosc for a particular problem. 





As we observed in Example 9-2, J, behaves fairly well, considering that parameters can assume 
only integer values. We conjecture that if 2, and n, were continuous, then a unique relative 
minimum for /) would exist. 
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| . 10. METHCD OF COMPUTATION OF THE TABLES i 

As indicated in the Foreword to this Handbook, the function tabulated, ; a 
a ey 
n | 
P(c.n.p) = > ( \px a (p+q=1) (10-1) _ 
x = ¢\X 
is related te the Incomplete Beta-Function Ratio | 


p 1 
I,m) = soo | Gas -! asi f sev (—sy"-' ds (10-2) 
ct) ct) 


by means of the relation 
Plen,p) = I Con eth) (10-3) 


On the other hand, the Incomplete Beta-Function Ratio Lfem), is known (Pearson! ) to satisfy 
the recursion formula 


1,(cun) =p Ise -— lm) tq Icom -— 1) (10-4) 
i:xpressed in terms of the function ? = P(c.2,p) this recursion becomes 
Plonp) =p Ple - In — Vp) +q Plen - Ip) (10-5) 


Eq. 10-5. together with 


P(O,n.p) =(p +q)" = 1 
and 





PQanyp) = p" 


permits the determination of all the tabular values desired. For fixed p, each value of P may be 
obtained for each n by increasing c from 0 to ” once all the values of P for n — 1 are known. Thus, 
for example, for p = 0.40, #2 = $2, since ' 











P(23,51,0.40) = 0.2723344, and 
P(24,51,0.40) = 0.1872942 are known, then ic 
P(24,52,0.40) = 0.4(0.2723344) + 0.6(0.1872942) = 0.22131028 F 
The computations for the original table, Army Ordnance Corps Pamphict ORDP 20-1, dated : 
September {952, were carried out on the IBM Relay Calculators. The values of P were computed to ie 
twelve decimals, which were rounded back to seven decimals for the final listings. For the x 
supplementary table--pages 578 through 626—involving p = .001(.001).G10, the Card Programmed 3 
Calculator, which was capable of operating with ten decimal numbers, was used. a 
XXXi : 
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